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ABSTRACT 


The use of boiling heat transfer in the liquid immersion 
cooling of electronic components has always been hampered by 
the excessive superheat necessary to initiate nucleation in 
the fluorinated hydrocarbons used as dielectric cooling 
fluids. In an attempt to overcome some of these difficulties, 
an experimental study of the effects of nucleate pooling 
boiling, on the heat transfer from surface near the boiling 
surface was conducted. An experimental chamber was constructed 
which had a column of four horizontal wires spaced 2.5 cm 
vertically from each other. The lowest wire was progressively 
heated from the natural convection region through the nucleate 
boiling region. A study was made of the effects of the boiling 
wake from the lowest wire on heat transfer from the upper 
wires. Under certain conditions heat transfer enhancements of 


up to 30% were obtained. 
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UNDTTsS DESCRIP LION 
m* Area 
Lak Thermal Expansion Coefficient 
J/kg Specific Heat 
m Diameter 
m/sec” Recelerationvol Gravity 
ND Grashof Number 
W/m*k Beak Traneter Coctricient 
W/mk Thermal Conductivity 
m Heater Length 
W Heat Transfer 
W/m* Heat Flux 

Prandtl Number 
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GREEK SYMBOLS 
Ns/m* Dynamic Viscosity 
m/s Kimemat ue Viscosity 
N/m Surface Tension 
kg/m° Density 
Various Uncertainty in the Variables 
SUBSCRIPTS 

proure 
Vapor 
Saturation 
Subcooling 
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I. INTRODUCTION 


A. BACKGROUND 

Research on pool boiling in electronic cooling systems has 
focused on three primary areas: (1) reducing the temperature 
excursion at incipient boiling, (2) reducing wall superheat 
during nucleate boiling and (3) enhancing critical heat flux 
(Mudawar, 1990). Since the heat transfer potential of nucleate 
boiling is well known, this regime was one natural selection 
for examination as a possible solution to this cooling 
problem. Nishikawa (1980) notes that there are two primary 
methods used to promote nucleate boiling. The first is to 
treat the surface in amanner that reduces its wettability. 
The second, more promising one, 1s to manufacture a surface 
with cavities which trap and hold vapor keeping the nucleation 
Site active. 

In recent years the drive toward higher heat flux and 
power density in electronic equipment has turned attention to 
the use of dielectric liquids. The 3M Corporation has 
commercially produced its "Flourinert" series of high 
arellectric, inert electronic cooling liquids. In contrast to 
water and other conventional coolants these liquids possess 
relatively low thermal conductivities and boiling points and 


small heats of vaporization. They demonstrate high wettability 


on most surfaces. Their surface tensions are the lowest ofivan 
known liquids, resulting in near-zero wetting angles. This is 
Significant because it eliminates many potential nucleation 
Sites. High dielectric strength and high resistivity are some 
of the characteristics that make these liquids well suited for 
electronic applications. Especially, FC-72 is known as 
chemically stable and leaves no residue during boiling. 

However, when using dielectric fluids, incipience 
superheat excursion 1S pronounced which leads to very high 
wall temperatures. A firm understanding of boiling incipience 
is required in order to implement effective cooling. Boiling 
can produce very large-scale heat transfer coefficients but it 
also introduces some problems: (i) Boiling restricts the 
physical design of the system. (ii) A high degree of superheat 
may be required if the surface is very smooth in order for 
nucleation to occur. Due to the high local heat flip 
occurrence of periodic boiling may result in local turbulence, 
even damage of electronic components. (111) Impuripres 
concentration may result at the nucleation site. 

Subcooling offers many advantages over saturated boiling 
for cooling large arrays of computer chips due to the 
decreased diameter and rapid collapse of departing bubbles. 

In this study heat transfer enhancement due to bubble 
pumping in FC-72 was investigated. It is hoped that they will 
serve to identify some trends and provide a basis for further 


research on actual components and integrated circuits. 


B. PREVIOUS WORK 

Hsu (1962) proposed a method for determining the range of 
active nucleation sites as a function of temperature or heat 
flux. This model revealed that the maximum and minimum sizes 
Semacelvye Cavities are functions of subcooling, system 
pressure, physical properties and the thickness of the 
superheated liquid layer. 

Bergles and Chu (1979) compared the nucleate pool boiling 
heat transfer characteristics of three copper Union Carbide 
Pol Fink test sections to a plain copper tube in distilled 
water and R-113. The experiments were conducted with three 
different treatments to the test surface prior to collecting 
gata: SubcoOoling; aging by preboiling in the pool; and heating 
the surface in air to remove all liquid in the pores. This 


research indicated: 


ie Heat transfer coefficients for the porous surfaces 
were four to ten times higher than for the plain 
Surfaces. 

Vig There was a significant temperature overshoot 


Prior to the initiation Of bolling and it was) oot 
very sensitive to aging or the power increment 
changes in water. These overshoots were sensitive 
to aging, initial subcooling and power increment 
Changes 1 R=11 35% 
Bergles et al. (1968), conducted nucleate pool boiling 
studies in water, R-113 and FC-78, using a 304 stainless steel 


tube. As heat fluxes were increased, the heat transfer 


coefficient followed the convective curve until a high 


Ssuperheat was attained, and then a distinct increase in the 
heat transfer coefficient was noted as nucleation occurred. 
This was followed by an abrupt drop in all temperatures by as 
much as 27°F. It was also noted that the incepritonmie 
nucleation could be triggered by vibrating the test surface. 
The conclusions drawn were: 

The phenomenon of temperature overshoot hysteresis in 


ordinary liquids is due to two causes: 


ils The existence of metastable bubbles which are 
triggered only at suffteient im high distureawe— 
levels, 

Lae The deactivation of larger cavities by displacing 


the vapor bys liquid duimmng subcooling= 


Many of the experimental reports published prior to 1986 
are reviewed by Bar-Cohen and Simon (1988). They noted 
Significant variations in the incipience Superheat from one 
report to the next. The authors discussed possible mechanisms 
for delayed nucleation. Recently non-repeatable and unsteady 
behavior of boiling incipience for a highly-wetting dielegeame 
liquid (R-113) was documented bysvYou vet ali. (1939)e 

The bubble behavior from the nucleate boiling regime to 
the film boiling regime has not been clarified yet, except for 
liquid helium. The transient boiling experiment of liquid 
nitrogen at stepwise heat generation in a platinum wire was 


carried out by Kunito Okuyama (1989). 


McAdams et ai. (1949), reported the effect of dissolved 
air on forced-convection subcooled nucleate boiling of water 
where the dissolved air reduced the heat flux at boiling 
PalGiacion. 

Torikai et al. (1970), used platinum wires immersed in 
water to measure dissolved gas effect on the incipience wall 
Supetneat under reduced pressure conditions (0.05 - 1.0, Bar). 

Surface enhancement is a very popular tool for improving 
pool boiling performance. Marto and Lepere (1982) found three 
commercially available enhanced boiling surfaces which shifted 
the boiling curve of FC-72 toward lower superheats compared to 
ee leainm Copper Surface, In- a study by Berenson. (1962), 
nucleate boiling heat transfer coefficients varied by 600% due 
to variation in surface finish. He emphasized the importance 
of the surface roughness. Recent investigations on surface 
effects in pool boiling by Chowdhury and Winterton (1985) show 
EMat Neat transfer appears to be unaffected by contact angle. 

ime “ciiecrs “or tlow velocity and “subcooling “on “eritical 
Meee Llux from a simulated electronic chip attached to the 
wall of a vertical rectangular channel was determined by 
Mudawar and Madox (1988). They concluded that heat flux was 
Significantly enhanced by increasing the degree of subcooling. 

ie Solubility of qas in hegqhnly-wettamng dielectric fluids 
moot Special concern: These fluids can.contain Up to 25 times 


more dissolved air by volume than can water: 48% air by volume 


can be dissolved in FC-72 compared to 1.9% air by volume in 
water (3M Product Manual, 1987). 

Experimental results presented by You, Simon and Cohen 
(1990) indicate that dissolved air plays an important role at 
incipience only when a very high dissolved air content is 
maintained. The surface tension value for FC-72 was estimated 
to have a 2% reduction due to high dissolved air content. 

S.M. You et al. (1989) conducted an experiment in 
Saturated R-113 fluid under atmospheric pressure with two 
types of wires. One is 0.13 mm diameter chromel wire (90% 
nickel) and the second is a thin platinum film on a quartz 
rod. They examined many pits of various sizes using a scanning 
electron microscope (SEM). 

Bar-Cohen (1990) conducted experiments providing pool 
boiling characteristics for a number of commercially available 
inert liquids (FC-43, FC-72, FC-75, FC-77, FC-84, FC-87) @auam. 
were obtained with a classical 0.25 mm platinum wire heating 
element. He concluded that: 

ala Nucleate boiling characteristics of these fluids 
depend on the heating surface, 


La Saturated pool boiling curves were almost 
identical with one another on a particular heater, 


S Superheat excursions at boiling incipience appear 
to always occur, but their extent can vary widely. 


C. OBJECTIVES 


The objective of this work was to investigate the effects 
of the boiling wake from a heated wire on the heat transfer 
from wires placed within the wake. This work was carried out 


Heing the dielectric liquid FC-72. 


II. EXPERIMENTAL DESIGN 


A. DESCRIPTION OF COMPONENTS 

Figure (1) shows the system configuration and identifies 
all major components. Figure (2) 1s a photograph (ore 
apparatus. 

The chamber consisted of a plexiglass box with an aluminum 
coverplate fitted with a rubber O-ring. This vessel has side 
wall heaters to allow the liquid to be preboiled for degassing 
before the run, and to maintain the liquid at Saturation 
temperature throughout the run. The vapor was condensed and 
returned to the chamber by gravity from the aluminum cover 
plate condenser equipped with thermoelectric coolers. The 
following is a description of the test apparatus. 

1. Test Chamber 

The test chamber was constructed of plexi-glass of 1/2 
inch thickness and had inside dimensions of 6 inches height, 
6 inches length and 2 inches width. An isometric view of the 
chamber is presented in Figure (3). Side walls are glued to 
each’ Obher- 

It has a 1 inch diameter hole and 1 inch protruding 
pipe and a copper elbow tube in the back face which is the 
pathway for the power and thermocouple wires. The chamber also 


allows the test section to slide in and out for the 


experiment. Cross-sectional and schematic drawings of the test 
section in the chamber are presented in Figures (4) and (5) 
respectively. The chamber also contained two small (1/4 in) 
pressure equalization outlets, located at the top-right and 
the bottom-left of the front vertical surface. It is a simply 
a closed tank with an aluminum cover plate. 

The test chamber is attached to a wood base in order 
to give better arrangement for cables and connections. 

2. Thermocouples 

Thermocouples, 0.01 in (0.254 mm) in diameter, were 
used to take measurements. The thermocouples for the plate 
condensers (thermoelectric coolers which are located on top of 
the aluminum plate) were bonded in grooves by placing a small 
drop of Eastman 910 adhesive on the thermocouple bead. The two 
thermocouples then were placed in the designated line groove. 
Using a weight, pressure was applied until the adhesive set. 
After allowing for three hours cure time, the remainder of the 
groove was filled with Omega Bond 101 Epoxy and smoothed to 
the level of the aluminum plate inside face. The Omega Bond 
101 Epoxy experienced little shrink after curing. 

3. Heaters 

The side wall heaters were glued on either side walls 
of the chamber. The actual size of the heaters were 125 x 11 
mm and they were subjected to 5 V, 2.5 amperes. The power load 


depended on the combination of the wires being heated. 


4. Heat Exchangers 
The three thermoelectric heat exchangers measure 1.56 
x 1.56 inches and were placed on top of the aluminum cover 
plate. The power supply for these were determined by the 
related curves. The property curves and technical data are 
presented in Figures (6) and (7) respectively. 
5. Aluminum Cover Plate 
An aluminum plate of 1/2 inch thickness was used as a 
cover plate for the test chamber. The inside face of the plate 
had two thermocouple beads to monitor the condenser 
temperature. Three thermoelectric condensers were placed on 
top of the aluminum plate. A thin groove was cut around the 
edges of the inside face for a O-ring gasket. The plate was 
screwed over the test chamber with eight stainless steel 
screws of 1/16 inch diameter and compressed the O-ring to form 
a watertight seal. Figure (8) shows the cover from the top 


with the condensers on it. 


B. INSTRUMENTATION 

Separate voltage taps across the heating element and 
across the precision resistor are used to compute heater 
resistance and heating element power. The heating element 
temperature is computed from its resistance via temperature 


calibration curve obtained prior to the experiment. 
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1. Power to the Heaters 

Each heater was run in series with a precision 
resistor that was measured to have a 0.1 ohm resistance. All 
four heaters were in parallel with a Hewlett Packard model 
member o739A, O-40 volts, O-1l.5 amperes, DC power supply. Both 
the voltage across the precision resistor and the chromel wire 
were measured independently. The current to each heater was 
calculated by the product of the heater voltage and heater 
current. Both voltages were measured by a Hewlett Packard 
model 3852S data acquisition system containing a Hewlett 
Packard model 44701A digital voltmeter, all controlled by a 
Hewlett Packard model 300 computer. 

2. Temperature Measurement 

The thermocouples for. the bulk fluid, ambient, ice 
reference point and coolers inside plate temperatures were 
referenced individually to an electronic Ice Point Reference 
as seen in Figure (9). Each reference thermocouple was 
connected such that its constantan lead was connected to a 
constantan lead of the measurement thermocouple. The copper 
leads of each measurement and reference thermocouple were 
connected directly to a Hewlett Packard model 44705A 20 
Channel relay multiplexer and inserted into the data 
acquisition system. The data acquisition system then measures 
the ice referenced voltage of each thermocouple. The voltages 


were converted directly into temperatures in the controlling 


a 


computer program by using a fourth order polynomial, fim 
the thermocouple manufacturer’s calibration data for copper- 
constantan thermocouples (Hazard, 1986). The locations of the 


thermocouples in the chamber are presented in Figure (10). 
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Photograph of Test Apparatus. 


Figure 2. 
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Figure 3. Isometric View of the Test Chamber. 
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Figure 5. Schematic Drawing of Test Section in Vessel. 
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Figure 6. Voltage Determination for Thermoelectric Coolers. 
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Figure 10. The Locations of the Thermocouples in the Chamber. 
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III. EXPERIMENTAL PROCEDURE 


A. NORMAL OPERATION 

In carrying out the experiments, the procedure was as 
follows: 

1. The ice point reference device was turned on. A mercury 
thermometer inserted into the ice base ensured that the ice 
Peeieeus at © + 0.05°C within “one: hour. 

2. The data acquisition system was then turned on. The 
internal voltmeter in the data acquisition system requires one 
hour warm-up period. 

3. The computer and printer were turned on. The data 
acquisition program contained-in Appendix D was loaded. 

4. The HP 6289A, 0-40 volt, DC voltage supply was set to 
zero and energized. The current limiter on the 40 volt power 
supply was set to the maximum current level. 

A sampling of all thermocouple temperatures was performed 
and printed. The power was increased incrementally and the 
output from the power calculation program was checked. 

The bulk fluid temperature as given by the wire resistance 
must match the value given by the bulk fluid thermocouples. 
When the results showed that the temperature stratification to 
be less than 0.10°C the experiment was ready to begin. The 


power to the test section is then increased by controlling the 
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test section heater voltage. The voltage was increased in 0.2 
volt increments until the initiation of boiling occurredyeae 
which time 0.1 volt increments were made until the maximum 
attainable heat flux was reached. This occurred at 5.7 volt 
for the power supply, making uSe of the HP 6289A, 0-40 volt, 
O0-1.5 Amperes DC voltage supply up to 8.1 volts, which 
corresponds to the maximum value. The power supply has a reset 
function to protect the device itself. The power is then 
decreased in the same fashion. At each power setting the 
system was permitted to stabilize for a few minutes and the 
following data were recorded: The heater voltage, precision 
resi Stor voltage, two thermocouple indi cateions Oe 
thermoelectric cooler temperatures, and three thermocouple 
readings of pool temperatures. This procedure was repeated for 


CVeLry "Seu, Ups 


B. DATA REDUCTION 

The test wire temperature was found using the calibration 
formula between thermal resistance and temperature value of 
the chromel wire. 

R=R,(1+e€AT) 

where, 

ANGE a d cae 2 

T, is the wall temperature 


T, 1s the bulk liquid temperature 
—€E is the slope of the calibration curve. 
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Details of the calibration slope value are presented in 
PepenOla sew hiss sendy mea lizes the walue ~or 0.0004 C/ohm 
offered by Kaye (1976). 

The starting value of the wire resistance is accepted as 
ene Ro in the previous formula. 

The heat flux is calculated from the following formula: 

Flux = Power/A 
where, A =n" DLis the total wetted surface area, and D and 
L are diameter and length of the wire respectively. 

The bulk temperature of the liquid was measured by taking 

the arithmetic average of the three thermocouple readings. 


es ote 8s) a 


The Rayleigh number based on diameter is given as: 


ODP Gea) eo 
Ra = eee 
VO 
where, 


Kr 
(p,-Cp) 


All the properties in the above equations are calculated 
at the bulk temperature. 


Peter nanatleon-of Nusselt number 1c: 


where h is the heat transfer coefficient and was determined 


ECOmM 


Zo 


h _ q” 


- se, 


where, 
T, = wire surface temperature (°C) 
h = heat transfer coefficient (W/m* - C) 
G = heat fuse /m) 


Churchill and Chu (1975) have recommended a single 
correlation formula for the natural convection from) )Yome 


horizontal cy anders. 


= ORS CWeR ani: 
Ming = (1012 6 0° 
[it (0. 55 oe a 


It is valid for a wide Rayleigh number range 
HO" Rae ese 


where the Prandtl number (Pr) is defined as 


AS 


IV. RESULTS AND DISCUSSION 


Data were taken and the related graphs can be grouped for 


discussion purposes as: 


ies individually powered wires 


Zs pumping effect runs 


A. INDIVIDUALLY POWERED 

Before investigating the pumping effect of the bubbles, 
individual boiling experiments on the wire heaters were 
@omaucted. 

Figure (11) shows wire 1 (the lowest wire) where distinct 
Single phase and nucleate boiling regimes are observed. When 
the heat flux was increased in small steps, the single phase 
Curve reaches 17°C superheat. At about 20°C an abrupt 
temperature decrease occurs as nucleate boiling begins. 

Figures (11) and (12) are essentially the same. The 
coordinates are heat flux versus surface temperature of the 
chromel wire. Figure (12) is on the log-log scale while Figure 
(11) is semilog graph. 

Figure (13) is the single run boiling curve of heater 
number 3. It was the only one to show overshoot of less that 


[os 


Va | 


Figure (14) for wire number 1 shows the same type of 
behavior as Figure (11). In this case data were taken at 
different times in order to demonstrate the repeatability of 
the experiments. Data points for two different runs did not 
coincide, but that is because the bulk fluid temperature is 
not the same for all experiments, but depends on the ambient 
temperature of the room. 

Figure (15) indicates all individually powered heaters on 
the same graph. The difference in the overshoot may be because 
of the surface geometry of heater number 3. 

The large aspect ratio of the heaters resulted in minimal 
conduction losses to the supports and yielded a uniform 
temperature along the length. 

The boiling curves for all four wires are practically 
Superimposable, even at the departure from nucleate boiling. 
They differ only in the heat flux at which boiling incipience 
occurred and at which full nucleate boiling began. 

Figure (16) shows the best fit to the natural convection 
portion of the experimental data for wire number 1, and is 
compared with Churchill (1975) correlation values. There is a 
constant offset from the correlation values because of the 
wall efffect except for the low Rayleigh numbers of that run 
which also shows considerable scatter in the data. 

The saturation temperature of the FC-72 is about 56. 
Teat ~ Ibnzxr 1S the subcooling. Figure (17) shows the eftfeces gem 


subcooling. At boiling incipience, a liquid mass near Mage 
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heated surface is slightly superheated, even though the bulk 
fluid is subcooled. Each run began with zero heating. The heat 
flux value was increased incrementally, continuing up the 
Single-phase natural convection line until boiling was 
initiated. The figure indicates that as the subcooling 
increases the curve shifts to the left, displaying significant 
increases in the overshoot and heat fluxes. With bulk liquid 
at 15°C, it shows large amount of Superheat, probably a bubble 
begins to form on the wire and as soon as the bubble grows a 
bit it encounters subcooled liquid and collapses. The liquid 
is so cold that it doesn’t allow the bubble to grow. As the 
bulk temperature of the liquid gets warmer, it allows the 
bubbles to grow. Subcooling results in smaller bubbles, which 
are less likely to spread to neighboring nucleation sites than 


the larger bubbles associated with saturated boiling. 


B. PUMPING EFFECT 

For this type of experiment, heater number 1 was used as 
the pumper. 

Wire number 4 was fixed at various heat flux values, while 
heater number 1 was increased in a stepwise manner. 

Figures (18) and (19) show the pumping effect of the 
bubbles generated by wire number 1 on the temperature of wire 
number 4. Clearly Figure (19) gives better understanding by 
way of marked corresponding points. As the temperature comes 


closer to the temperature at which the boiling incipience 


Zo 


Starts, wire number 4 is little effected. Just after wire 
number 1 is subjected to overshoot it also has similar a 
experience, i.e., a sharp jump . Actually, it was the general 
pattern for most of the other runs. This Suggests that boiling 
can be started by the pumping effects of the rising bubbles 
which trigger the heater located above it. The characteristics 
of the bubbles leaving the lower heater are uncertain. They 
may be the wake of the buoyancy-induced motion, especially 
when the upper wire was fixed at a high heat value. A similar 
trend can be observed from Figure (20). One thing to notice is 
that heater number 4 is fixed at a heat flux that corresponds 
to a starting value of surface temperature of 50°C for this 
run. This is below the temperature that meets the overshoot. 
Less jump in the temperature is associated with the lower 
fixed heat flux value. If we compare it with Figure (19), 
which was at higher heat value than this one, Figure (21) 
provides a detailed look at Figure (20). 

Figure (22) shows wire number 3 fixed at 37,490 W/m aeeee 
has the points marked in chronological order. The 3rd wire 
heater is located closer to the first wire than is the fourth 
wire, so it allows most of the bubbles to reach the upper 
heater with collapsing due to surrounding subcooled liquid FC- 
Ze 

Figure (23) shows another combination of heater number 3 
with heater number 1. Heater number 3 is kept at 20,450 W/m’. 


This was very low heat flux value that the triggering effect 
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Gid not show itself. Bubbles still caused some decrease in the 
temperature of the upper wire. 

Figures (24a) and (24b) show very similar results. In this 
case wire number 2 is kept at 20,530 W/m’. Since wire number 
3 is the closest to wire number 1. Being so close helped it to 
lose some heat corresponding to a decrease in its temperature. 

Figures (25a) and (25b) show no effect from the bubbles 
from below. Here heater number 2 is kept at a very high heat 
mero ot 2i2, 700 W/m, which is already boiling. It is not 
sensitive any more to the bubbles produced by heater number 1. 
Consequently no jumps or apparent decrease on the surface 
temperature of the number 4 wire are seen. 

Figure (26) shows similar trend as the early graphs. In 
this case the bubbles are produced by wire number 2 and the 
effects on wire number 3 are presented. Heater number 3 had 
the temperature decrease as the pumper had overshoot 


experience. 
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Figure 11. Wire Number 1 Individually Powered, Semi-log. 
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Figure 12. Number 1 Heat Flux vs. 
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Figure 13. Number 3 Shows Less Overshoot. 
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Figure 14. Number 1 at Different Times. 
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10° 
an 


@ 

e si 

= A & 

oe a gO 

aS 

3 — Se ae ee 
> Reo 

(23 aaa Sa 


LEGEND 
a = WIRE #1 
o = WIRE #2 
4 = WIRE #3 


+= WIRE #4 





2020 34.0 45.0 56.0 67.0 78.0 89.0 
TSURFACE (C) 


TBULK=23.5 C TSAT=56 C 
Figure 15. All Wires Individually Powered. 
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Figure 16. Best Fit Experimental Data Compared with 
Churchill (1975) Correlation. 
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Figure 17. Subcooling. 
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10° 


i 


HEAT FLUX (W/M 2) 


LEGEND 
1 INCREASING 





23.0 34.0 45.0 96.0 67.0 78.0 89.0 


TSURFACE (C) 
21JAN14G 


Figure 18. Pumping Effect with Wire Number 4 at 42,800 W/m’. 
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Figure 19. Pumping Effect with Wire Number 4 
(Detailed Figure). 
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WIRE #4 AT CONST. HEAT FLUX 
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Figure 20. Wire Number 4 at Constant Heat Flux (38,600 W/m’). 
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Figure 21. Detail of Figure 20. 
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WIRE 1 & 3 (FIXED) 
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Figure 22. Wire Number 3 at Constant Heat Flux Value of 
37,490 W/m’. 
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Figure 23. Wire Number 3 at 20,430 W/m’. 
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WIRE #¢ AT CONS. HEAT FLUX 
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Figure 24a. Wire Number 2 at 20,530 W/m’. 


45 





WIRE #2 AT CONS. HEAT FLUX 





a ADOROW wh te O 


ae 


(2 WAN) XO'D BVO 


TSURFACE (C) 


Figure 24b. Detail of Figure 24a. 
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Figure 25a. Wire Number 2 Kept at 212,700 W/m’. 
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Figure 25b. Detail of Figure 25a. 
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Figure 26. Wire Number 2 with Wire Number S7ac 
Constant Heat Flux. 
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V. CONCLUSIONS 


The figures give the idea of the pumping effects. When the 
pumping heater is powered progressively, the upper heater will 
react as indicated in the following. 

° If the upper wire is located close to the lower 
pumper, it will meet most of the bubbles. It will lose some 
heat depending on how close it is. 

° The initial heat flux value of the upper heater is 
also important. The amount of temperature jump which 
Simultaneously occurs when the pumper has overshoot depends on 
how warm the upper wire is. 

° If the upper heater already has bubbles, this means it 
had incipience before, and the bubbles will not help it lose 


more energy. 
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VI. RECOMMENDATIONS 


1. Use platinum wires rather than chromel wires to get 
better temperature measurement. 

2. Before the start of each run, boil the pool for a 
sufficient time to degas the liquid. 

3. Modify the software for taking data by adding a WAIT 
Gommand to obtain the time function. 

4. An addition of a foil type heater placed horizontally 
at the bottom of the chamber may help create a uniform 
temperature throughout the bulk fluid. 

5. An extension of tygon tubing of about one foot length 
would take care of overflow and hydrostatic head determined by 
the liquid level inside the filling tank. 

6. The present study deals with only subcooling effects. 
Saturated pool boiling should also be studied. In order to 
have constant temperature in the bulk fluid, the thermo- 
electric condensers and the bulk heaters can be run 


Simultaneously. 
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APPENDIX A. SAMPLE CALCULATIONS 


A. DETERMINATION OF INPUT POWER 
For the sample calculations, the output values of wire 


number 1 were used: 


ViOws Game =) 
Omel 


where volt (12) and volt (11) are the voltages across the 


Power(1)= VOLe« lee) 


precision resistor and wire respectively. 


-( Bagee® 
Lo 


|(1.988)=0.96 watts 


B. CONVERSION FROM VOLTAGES TO TEMPERATURES (°C) FOR 

THERMOCOUPLES 

The HP 3852 A Data Acquisition System channels 1-9 record 
the voltage readings in millivolts. The following polynomial 
formula is for copper constantan thermocouples. 

T(°C) =0.10086091 + (25727.9 * emf) - (767343.5 * emf? ) 


+ (7802.5596 * emf?) - (9247486589.6 * emf’) + 
(6.98F11* emf?) - (2.66F13 * emf®) + (3.94E14 * emf’ ) 


Using thermocouple 1 


emf = 0.000779 volt To) 2c 


C. CALCULATION OF WIRE TEMPERATURE 
The calibration curve for the chromel wire (Kaye, 1976) 


yields the slope € of 4.107 ohm/°C. 
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R=R,[(1+eET-T,] 


Here R, and R are the resistances of the wire at the reference 
point and at any point to be measured R, and ohm(1) are the 
corresponding variables used in the software. 


Mii = «(fonnii) 7k) =91) \/ Kaye + £ 


where, 
del) = wire temperature 
Kaye = 0.0004 ohm/°C constant slope 
Te = wire reference temperature (°C). 


T.(1) = (((4.1009/4.075) - 1) / 4E-4) + 19.64 


= S20 6 


D. CALCULATION OF THE AVERAGE BULK TEMPERATURE 

The average temperature of the ligquid was simply the 
arithmetic average of the three thermocouple values immersed 
mmoteme biti £luid. | 

Oaik = wn! ea (2) eo) ae 


SOF 6 954 a1 06) Se 1 9b ee 


E. DETERMINATION OF HEAT FLUX 


PEO s(t y= eee 
AREA 
where, 

AREA = ae a eo a (oe by ke 

De i = Diameter and length of the wire. 

AREA (5 1 415970) 1000S eon os als) = S2042-5 m 

HPLUX = 22222 = 77460 W/ m2 
-04E-5 


Sie 


F. FLUID PROPERTIES DETERMINATION 
1. Thermal Conductivity k(W/m°C) 
From Figure 5 of the 3M Corporation Flourinert Produeas 
Manual, the thermal conductivity coefficient curves have been 


determined to be 


W 
(0.6033 —0.00115 +m. ano aR 





at f, = 50°C it yields 02054587 m- 
2. Liquid Density p (kg/m’) 
Using the expression in Table 4B and constants 
presented in Table 4C of the Product Manual yields 
Seal 740 == 70 Oo) LOU 
T, film temperature must be in units of “CC. Atel, = 50°@eeame 
yields™p— lc0oesexc/m - 
3. Kinematic Viscosity v (m’/s) 
From Figure (3) and determining a 4th order curve fit 
yields 
v = (0.54 - 0.0116 * T, + 0.0002085 * T,’) * 1E-6 
at T, = 50°C this yields) Va——4 creme oe 
4. Density of Vapor p, (kg/m’) 
From Table 13A of the Product Manual, values for vapor 
density can be taken. 
For the range of 20-80°C the data can be approximated 
as a fourth order polynomial. 
Vi-= (246° —= 10. Ose ee oe ee 


0.00193 #19 F-61558 6a) ey oO 
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T, = 50°C it yields 


Oo, = Ge67 (coy ml: 


5. Specific Heat C, (J/kg°C) 
For all Pleurinere Elect rochemicals (Figure 4) 
Coe Ones a ae) SEs Te a ep 
emete— 50 CC; Calculation yields 1086.7 J/ka Cc: 
6. Thermal Diffusivity a (m’/s) 
a = k/PC, 
at T = 50°C, calculation yields 3.115E-8 m’/s. 
7. Thermal Expansion Coefficient $ (1/°C) 
Table 4B in the Product Manual shows 
Be 205002617 (12740. =.0 700261) 34 -7,) 
(wee oC. Walculation results ino 0.00162 .41/°¢C. 
8. Prandtl Number (Pr) 
Pr = v/O 
ceet=—o) Colhis yields Pr-= 15.4 
9. Calculation of Latent Heat hfg (J/kg) 
Table 13A of the Product Manual shows the values can 
be approximated by a linear equation. 
eG c= Oo. = ee 3 ok oe) ee OS 79 
mere 10,79 1s the conversion factor from Cal/Mole to (J/kg) ; 


At T = 50°C calculation yields 74672 (J/kg). 
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APPENDIX B. UNCERTAINTY ANALYSIS 


A. UNCERTAINTY IN AREA 


OQ, = .0001 inches 
QO 7=) += 020 sinenes 
A = 6 0) 
where: 
D = Wire diameter 
L = Length of heater wire 
A = (3.14159 35 4 502005 aes) 
= OP Oe lion. 
= | Stee On om 
OA 7 
Ob ae 
0A 7 
ab mie 
Therefore, 
®, = [(mL@,)? + (mDw,)?]°” 


®, a 2 2 
A 7 we + Me 
D re 


315 


1/2 
| sea — 
FOR oss SS 
on = 0.021 or 2.1% 
A 


a) SG 2 ea ee) Oe 


UNCERTAINTY IN POWER 


0), = +0.05 Ampere Taken from (Hazard, 
calculation) 

WO, = 20> VoLt 

OQ = Vi 

00 = 

vee ue 

dQ = 

Gi) °° © 

W, ae 1 ee Gi) eee Vie aye 2 
i ar 

ba 2 | 0 oS} =| 0 Ea 

O 7.46 0.018 


OOS Oi war > 


a 


1986) 


UNCERTAINTY IN HEAT FLUX 


le 


Q 
A 


i 


[ 
Or 


02021 + 0 nts ss 


O62 5or see Ze 


aye 


iRroe 


D. UNCERTAINTY IN TEMPERATURE 


Using thermocouples: 


Variable Uncertainty 
Voltmeter resolution a0 26 7c 


1.0 mv resolution 


Ice reference bath GE05-C 
Temperature 
Polynomial Temperature O20CGG3°C- RMS 
Conversion 
R Iba. 

OT, = PRG ene CO eins eee 


ne O02 5) e005) aie) Oe 00 668 


QT, OU CS 1026 rS 6 


E. UNCERTAINTY OF WIRE SURFACE TEMPERATURE 


Prom the calibration formula 


i= Ra lee Ae | (1) 
where, 
AP =e wearer a7 
E = Slope of the calibration curve 
IRE = Bulk fluid temperature. 
Brom 1) 
dye = [eo Jeon POm els 44 Ace 


Sy, 








sar . = 2em0'G 

Li b 

ee = 1 © = 0.0004 + 12% 
= R.€ 

OT -R 

IE a R. = 3.88 Q 

i Roe < 


5 
ce 


€ 
La | 
] 


s { 1 (0:063)* 4+ (6.49)" 26627" = 


hy toil 


F. UNCERTAINTY IN AT 


Where AT = T, - T, 


60 


OAT 
Es 


War 


a a4 








Oar Z lt aor a7 -S 
OR 21 544.5), 4 — 23-99) 


Ao 


UNCERTAINTY IN HEAT TRANSFER COEFFICIENT 


q 
AT 


AT 


61 


q is the heat flux 


G= 2010 9CnW7 m- 


dh 
OAT 


c 


it 


lI 


AT = 44,14. 
2 5 LFZ 
ee | coe = 20.15°C 
AT g AT? ek 


Leyeg 
(0.062)? +( 9a ) 


US 22200 (22.6 


GZ 
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APPENDIX C. CALIBRATION OF CHROMEL WIRES 


The chromel wires of about 3 inch length were immersed in 
a Rosemount constant temperature distilled water bath, and 
were compared with a platinum resistance thermometer as a 
Standard. The temperature was raised from 30°C to 90°C in 
increments of about five degrees. At each increment, the 
temperature was allowed to stabilize for five minutes, then 
measured uSing a Newport digital thermistor devise for the 
0.005 inch diameter chromel wire. All recorded data are 
presented in the following table. The corresponding best fit 
line is also presented on the graph. Two sets of data were 


jaken at different times. 


Thermistor T/C Resistance 
a anes OS aire ned AE 2 ke ae 

29 BObyZ 4.74 
Cy oer Sos 4.76 
49.5 ah 4.78 
61.9 63.4 4.8 

UoeZ ie oS 
OM] a 88.4 4.85 
D283 Bo 4.79 
Seo bo: 4.8 

60.5 62 4.8 

68.3 69.6 4.81 
Ou 72 4.82 
Ts oe Aa oi 
86.2 86. 4.84 
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Figure Cl. Chromel Wire Calibration. 
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APPENDIX D. SOFTWARE 


Software written in the BASIC language converts voltage 
measurements into temperature values for the bulk fluid, ice 
reference, and coolers inside the thermocouples. It uses the 
a value of 0.0004 Ohm/°C to calculate the surface temperature 
offered by Kaye (1976). The printed output covers bulk fluid, 
wire surface temperatures and heat flux. 

A FORTRAN program reads the output values created by the 
previous program and calculates liquid properties R,, Nu, and 
CoA. 

The GRAPHER program graphs the data and prints it uSing a 
he —-7440A Plotter. 

Initial experiment data was fed into EASYPLOT software on 
the Naval Postgraduate School’s main frame computer via floppy 


diskettes. 
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fu 
uf 
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SI wn 
Ol ha = 
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U 
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mn 
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oO 


OPTION BASE } 
BEEP 
PRINT "INPUT GATE” 
INPUT Tarihs 
PRINT " INFUT TIME " 
INPUT Zam 
'FRINT “SCANNING MEAS" 
REAL Emf(7),7(7) ,0¢9) Volt(tt:19) Power(4) Hflux(4),Ts¢4) ,Ohm( 4) 
REAL Sumo 
OUTPUT 705; "CGNFHEAS OGY ,1G!-107,USE @" 
FOR I=1 TO 7 
ENTER 709;:Emf( 1) 
iPRINT IsEmfi i) 
NEXT I 
OUTPUT 789; "CONFHEAS OCV ,11i-119,USE 2’ 
FOR I=11 70 19 
EHTER 7OSsVoltii i 
(PRINT isVoit( i> 
NEXT I 


Bi i d=. 108 =e 


FOR J=i TO & 
SumosSumoto0( J )etmF( la (d-i 2 
NEAT J 
T(1I3=Sumo 
Tost 1(Z IFT S)t11 1 75 
NEXT I 
PRINT "GATE + +larins 
PRINT "TIHE 3; "sZam 
PLN oie Bet tise UPPER PORTIGN GF THE FLUID” 
[PRINT USING “1°43 80 (bert ye ee UPPER PORTION OF THE FLUID" 
PRIN) 0 ogee Ieee LOWER PORTION OF THE FLUID” 
PRING “T3529 s 7s" oC LOWER FORTION OF THe FLUID” 
PRINT Tosa Cad ot TE COGLER INSIDE AL. PLATE” 
C 
C 


TT) 


esc) US Tie TR Oy ra TE CUULER INSiBe ALE enie 
PAN le Wt oor brrs AMBIENT (ROOM TEMPERATURE )" 
TENPERROURE 7° 
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wei a) @w@ =3 GC) GO) Gi oD 


“POPS Om - @ @ OO & Or 


oO «oO UW) UW) 


GS & w wm 
— 6) 


iz 


BRN nl 37.1 t pu cnaG ICE (REFERENCE TEMPERATURE)" 
INEXT I 

PRINT "“TEULK="3Tb 

Ipower coiculationsr+ssssessess 

Fower( |) =(Volt(iZ)/. 1 deVoltC Try 

Power (3 9=(Volt(17)/.1)*Voitt ib) 

theater surface area calculations 


LRRRETERAHREHREEREE EEO ERO PERE EDRESD 


% 


Bia=.805*.0254 
Aiw3*.0254 
Alan=(3.1415526 }*DiarAl 
lheat flu» cale. 
LPH HSH RROERROTOARROEOR®E 
Hriux(t =Power( 1 )/Alan 
Rrlux 2 )=Fower(2)/Aian 
Hfiux(3 d=Powert( 3 i/Alan 
PRINT “powerl=";Fower( 13 
IPRINT “powerZ=";Fower({ 2) 
IPRINT "power3=";Power(3} 
[PRS TRS O HEFT SHAPER RAS EREASD 
GnmCh p=Volt( li d/(Voitii2)/.1) 
IGhmi 2 =Voitlida/(Volttis)/.15 
GnmeSd=Voltc 1b d/iVoitii7s.132 

'PRINT "“OHMZ=";O0nmiZ) 
fPRINT "GHMi=";:Grmij) 
[PRINT "GHMS="s0hm(3) 
Ot=To-1a 
1 Ya iS THE TEMP OF BULK FLUIDS Tb WKEN FOWER OFF 
! Ro IS THE RESIS. OF CHROMEL WIRE ACTUALLY I7 is 
i TRE Ohm( 1) VALUE WHEN THE SULK FLUID TEMP (Tp=Tamb) AND POWER OFF 
| ( POWER OFF MEANS Power(1?76.0084 and it is obtained by setting ) 
i POWER(4)°O.G088 
1 RQl=4d.069 RO4=3.7981 

(Vic 40 V Supplyist to 30.1 U.i0 position WHEN 1&4 INSERTED 

Kaye is the value of the calibration slope Refis.M YOU et al} 


’ 
{ 
i SHES THA H ETH SERRATE RE HEEE RE ERETEREEHREEREEERRRE DERG 
= 
j 


a=TD 
Roi=z4.9065 

Rofs=5.7S5i 

Kaye=,0004 

R=4.959*(1+.0004-D0t 3 

Ts¢ i o=d¢ COnme i /Roi )-1 )/Kaye t+Te 
TstS oC ( (Gh 3 9/Rod 1 b/Kaye tia 


1335 Trilm=(Tb+Ts(4))/2 

roe PRINT “TS i=" sis 7° RPLUS 73 Ht luxe b) 
is RIND eoTSss=" ists)” HELUG wa star xe oO) 
iS6 TRIRINT "T film s"s7fiim 

340 END 


Ge) 


** * EK € 


¥ 


4b 4b 646 3 eb i Fe ET FET TA HE TEE IE EEE TE EI He Ee ttt 0-10-1018 -10-28-18-40-4--40- 9p 
PARAMETER ( Nees ) 
INTEGER I 
REAL TS (N), TB(N) , HF (CN), NU(N), RACN), NUC (N), TE (N), TE1(N), BETA(N) 
+ » DENF (N) , DENY (N),KINV(N),CP(N),K(N) ,HFG(N) ,ALF(N), PR(N) 
rn » SURT (N) , RAYD (N) , RAYFB(N) , NUSEX (N) , NUSC (N), DIA, G, V(N) 
+ »MUC(N),Y(N),X(N),WS(N), TSAT, CSF (N) 
sheet tiehietetiatiahialehathetet ReteRehet tetet Mhetataitehm ttt ttt tte ce Oe eee 
READ Tsurface,Heat Flux, Tbulk FROM DATA FILE * 
WRITE Tsurface, Tbulk, Heat Flux, Nusselt, Rayleigh, Nusselt (Chur) * 
Tfilm , beta, density, viscosity,Cp, conductivity, latent heat, 
Prantl, surface tension 
SFE AE OG IE EE EE EE EE IE EE Ee EE EE 6-8 bb 1 ee 4 tb 3 3 6 3b 8 4 28 46-8 8 
OPEN (UNIT=1, FILE=’ INPUT. DAT’ , STATUS=" OLD? ) 
OPEN (UNI T=e, FILE=’ QUTPUT. DAT’ , STATUS=" NEW? ) 
OPEN (UNIT=2, FILE=* PROPE1. DAT’ , STATUS=" NEW? ) 
OPEN (UNIT=4, FILE=’ PROPE®@. DAT’ , STATUS=’ NEW? ) 
OPEN (UNIT=S, FILE=" ROHSE. DAT’ , STATUS=’ NEW’? ) 


DO 15 I = 1,N 
READ (1,*) TS(I),HF (I), TBCI) 

TF (I) =(TS(1)+TB(I)) 7/2. 

TF1 (1) =TF(I) +275.15 

BETA( 1) =O. 88261/ (1. 740-8. 8B261*TB(I)) 

PRINT*,’BETA IS =’, BETACI) 

DIA=B. 2O5#2. 9254 

G=9.81 

TSAT=56. 

WS (I) =TS(1)-TSAT 

DENF (I) =(1. 740-0. 89261*TB( I) ) #1202. 


V(I) = (246-18. 9S*¥ TBI) #2. G1 S# TB (I) *##2-19. SE-4*TE (1) 43+ 
+ 6.55E-6*TB(I) #4) /SS8. 

DENV (I) =1./V (1) 

KINV (I) =(2. 54-0. 8116*TB( I) +0. BOVCO8S* THC I) #¥2) #1 E-6 

CPCI =(0. 24111145. 7OS7E-4G#TE(1I) ) #4186. 

KCI) =(—-0. VOL1IS*TF (I) +8. 6055) /1. 

HFG (I) =(8057~-22. 7S*TB(I) ) #18. 79 


ALF (I) =K(I)/ (DENF (1) *CPCI)) 
PRCI)=KINV(I) /ALF (I) 

SURT (I) =(-8. 12*TF (I) +15) /1980. 
MUCT) =KINV (I) #DENV (TI) 


| BSE SG HE IR Lb ie ie ER Sah ih HE BR le Ub BRS iH He HE iS IR EE bb i bE ES a 
Y CI) =CHF (I) / (MUCI) HFG C(I) )) *SQRT(SURT (I) /(G*(DENF (I) ~-DENV(I)))) 
XC(QTISCP CIV #¥ (TS CII -TBCI)) /CHFG CI) ORCI) #1. 7) 
CSF CID =X CID /S (CV CI) ##0. S32) 

Se He Yee he ee Ie FR FE te IE i a HE FE Je JE i J I tb i tb i i FE IE Hb J Hb ek ib I i Sb J J J Ed te PE et ee 
CALCULATE RAYLEIGH and NUSSELT NUMBERS 
8 ee ee ee Dt te ee lle 
RAYD (1) =G¥BETA(I) #¥(TS(I) -TB(I)) #DIAw#S/ (KINY (I) ALF (ID) 
RAYFB (1) =G¥BETA(I) HF (I) ¥DIR#*#4/ (KCI) #KINY (I) HALF (I)? 
NUSEX (1) =HF (I) #DIA/ (KCI) * (TS (I) -TBC(I))) 
NUSC (1) =(@. 60+ (Q. SA7*#RAYD (1) #61. /6. 9/6 (1418. SSI/PR(I) ) HH C 

+ 9./16.)) we (8. /27.))) ) #2 


WRITE (2,20) TS(I),TB(I),HF (I), NUSEX (I), RAYD(I), NUSC(I) 
WRITE (3,20) TFE(I), BETACI), DENF (I), DENV(I),KINV(I),CPC(T) 
WRITE (4,40) K(I),HFG(I),ALF(I),PR(I),SURT(I) 
WRITE (5,50) Y(I),X(I),CSFC(I) 

FORMAT (1X,F8.2, 2X, F9. 1, 2X, F12. 2, 2X, FS. 3, 2X, F12- 7, 2X, F15.3 ) 


32 
40 
~9 
te 


FORMAT (1X,F8. 4, 2X, F9. 4, 2X, F9. 1, 2X, F8. 4, 2X, F18.5, 2X, FiS5.5 ) 
FORMAT (1X,F8. 4, 2X,F14. 4, 2X, E18. 3, 2X, F14. 4, 2X, F9.S ) 
FORMAT (1X,F8.4,2X,F14.4,2X,E18.3 ) 
CONTINUE 
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